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ABSTRACT 


Molecular  dynamics  simulations  were  performed  to  characterize  the  lattice  morphology  in 
the  region  adjacent  to  the  interfaces  in  nanocomposite  systems  of  a NijAl  dispersoid  embedded 
in  A1  matrix  (NijAl/Al)  and  an  AljNb  dispersoid  embedded  in  aluminum  matrix  (AljNb/Al).  A 
nearly  perfect  coherent  interface  is  obtained  in  the  ALNb/Al  system  with  the  lattice  planes  of 
dispersoid  and  matrix  aligned  parallel  in  all  directions.  The  simulation  results  show  the  presence 
of  the  matrix  atom-depleted  regions  near  the  dispersoid  boundary  for  most  cases.  Detailed 
analysis  revealed  that  certain  sites  immediately  next  to  the  dispersoid  are  energetically  favored 
for  the  matrix  atoms  to  occupy.  The  matrix  atoms  occupying  these  sites  attract  other  atoms 
producing  the  depleted  regions.  In  certain  specific  situations  of  AljNb/AI  system,  however,  the 
wetting  of  a rotated  dispersoid  is  overwhelmingly  complete  prompting  the  need  of  further  study 
for  better  understanding.  The  order  parameters  of  dispersoids  calculated  for  N13AI  in  aluminum 
is  nearly  constant  while  that  for  AIjNb  in  aluminum  is  rapidly  decreasing  function  of  temperature 
in  the  range  of  300  to  1800K. 

INTRODUCTION 

Second  phase  particles  that  are  mixed  and  uniformly  distributed  within  a matrix  are  called 
dispersoids.  By  selecting  appropriate  type,  size,  morphology,  and  concentration  of  these 
dispersoids  in  a matrix,  it  may  be  possible  to  develop  structural  components  that  are  capable  of 
actively  indicating  their  physical/chemical  states  and  are  optimized  for  best  performance  based 
on  given  application  purpose.  For  example,  the  magnetic  properties  of  NijAl  particles  can  be 
affected  by  the  local  stress  states  that  may  be  indicative  of  the  changes  in  other  material  states. 

At  the  same  time,  when  mixed  with  a substantial  volume  fraction,  the  small  nearest  neighbor 
distance  between  such  densely  populated  dispersoids  may  prevent  the  dislocation  loop  formation 
while  being  highly  resistant  to  cutting  by  dislocations  [1],  The  current  series  of  investigation 
concentrates  on  the  morphology  of  such  composites  to  establish  a foundation  for  systematic 
study  on  nanocomposite  systems  of  damage  indicating  capabilities  and  prolonged  service  lives. 

Our  previous  study  clearly  showed  that  a NijAl  dispersoid,  which  is  coherent  with  the 
aluminum  matrix,  causes  a substantial  degree  of  local  lattice  distortion  in  the  vicinity  of 
matrix/dispersoid  interface  [2]  that  is  known  to  enhance  the  strength  of  alloys  [1],  An 
unresolved  problem  of  the  study  was  the  appearance  of  lattice  atom-depleted  regions  near  the 
interface  in  the  composite  system  with  a rotated  dispersoid  and  its  origin  was  not  well 
understood.  Therefore,  the  purpose  of  the  present  study  is  to  find  ways  of  minimizing  the 
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regions  of  matrix  atom  depletion  near  the  interface  and  closely  investigate  the  nature  of  these 
regions. 

The  lattice  structure  of  Ni3Al  is  known  as  Llj  which  possesses  a cubic  symmetry.  Being 
paramagnetic  in  nanosizc  [3],  the  magnetic  field  effects  of  Ni3Al  particles  are  not  considered  in 
the  present  simulation.  The  stable  phase  of  Al3Nb  is  known  to  be  DO22  structure  which  is 
tetragonal  instead  of  cubic  Ll2  [4].  Nevertheless,  we  choose  the  mctastable  L 1 2 phase  of  Al-,Nb 
in  the  present  study  for  two  reasons.  First,  maintaining  the  cubic  symmetry  of  the  system,  the 
simulation  results  can  be  directly  compared  between  the  two  systems  of  Ni3Al/Al  and  Al3Nb/Al. 
Second,  by  varying  the  physical  and  chemical  states  of  the  metastable  system,  some  unexpected 
but  interesting  facts  may  be  found. 

DETAILS  OF  SIMULATIONS 


For  each  system  of  Ni3Al/AI  and  Al3Nb/Al,  a dispersoid  is  constructed  at  the  center  of  a 
cubic  block  of  pure  aluminum.  Simulations  arc  first  performed  by  arranging  the  lattice  planes  of 
dispersoid  and  matrix  parallel  in  all  three  axes  of  Cartesian  coordinates.  Then  the  dispersoid 
lattice  is  rotated  with  respect  to  that  of  the  matrix.  As  an  initial  investigation  of  a systematic 
study  on  nanocomposites,  the  rotation  of  the  dispersoids  in  the  present  study  is  limited  to  the  x-y 
plane  and  three  angles  of  rotation,  i.  e.,  15°,  30°  and  45°.  Also  in  this  study,  the  dispersoid  lattice 
spacings  are  initially  set  to  be  very  close  to  their  real  values,  i.  e.,  0.357  nm  for  Ni3Al  and  0.384 
nm  for  Al3Nb  compared  to  0.405  nm  for  pure  aluminum.  The  inner  radius  of  matrix  and  the 
outer  radius  of  dispersoid  are  chosen  to  allow  partial  overlap  between  two  lattices.  The 
overlapping  lattice  points  are  removed  to  ensure  a minimum  separation  among  the  atoms  in  the 
interface  region.  With  further  adjustment  of  these  radii,  the  two-dimensional  dispersoid 
boundaries  in  the  x-y  planes  arc  selected  to  be  Ni-only  edges  for  Ni3Al/Al  and  Al-only  for 
Al3Nb.  The  simulations  are  performed  in  the  temperature  range  of  300  -1800K  based  on  the 
temperature-independent  parameter  values  for  a modified  embedded  atom  method  (EAM) 
potential  given  by  Baskes  [5], 

RESULTS  AND  DISCUSSION 


Fig.  1 shows  the  results  of  simulations  performed  at  300  K on  the  Ni3Al/Al  system  with  their 
lattice  planes  arranged  parallel  in  all  three  directions.  Shown  are  the  atoms  in  three  planes  of 
(001)  centered  at  z = 0.  The  hexagons  and  tiny  dots  in  this  figure  represent  Ni  and  A1  atoms, 
respectively.  Similar  results  obtained  at  900  K arc  shown  in  Fig.  2.  The  presence  of  the  matrix 
atom-depleted  regions  around  the  dispersoid  is  clear  in  the  results  of  these  two  figures.  One  of 
the  goals  of  the  present  study  is  to  find  the  origin  of  such  matrix  atom-depleted  regions  at  the 
interface,  and  a simple  model  has  been  established  and  applied  to  the  results  of  Figs.  1 and  2,  as 
explained  in  the  following. 

The  dispersoids  shown  in  Figs.  1 and  2 are  almost  octagons  and  two  types  of  edges  are  found 
in  their  boundaries;  those  consisting  of  four  nickel  atoms  and  the  others  consisting  of  five  nickel 
atoms,  which  we  call  the  four-atom  edges  and  five-atom  edges,  respectively,  for  brevity.  Upon  a 
close  inspection,  one  finds  that  three  matrix  aluminum  atoms  are  lined  up  immediately  next  to 
the  four-atom  edges. 
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Figure  1.  Simulation  results  of  the 
NijAl/Al  composite  system  at  300K 
without  rotating  the  dispersoid. 


Figure  2.  Simulation  results  of  the 
NijAl/Al  composite  system  at  900K 
without  rotating  the  dispersoid. 


Fig.  3 shows  lattice  sites  near  the  interface  for  the  matrix  aluminum  atoms  to  occupy. 
Referring  the  results  of  Figs.  1 and  2,  it  is  clear  that  the  sites  A and  A’s  next  to  a four-atom  edge 
are  energetically  more  favored  than  B and  B’  which  are  next  to  a five-atom  edge.  In  Fig.  3,  one 
can  see  that  an  aluminum  atom  in  site  A has  four  dispersoid  Ni  atoms,  denoted  by  x,  at  the 
nearest  neighbor  distance.  The  same  applies  to  an  aluminum  atom  occupying  site  B.  The 
important  difference  is  that  a complete  two-dimensional  unit  cell  is  formed  when  an  aluminum 
atom  occupies  site  A.  Once  this  occurs,  A’  is  the  next  favorite  site  for  occupation  since  there  are 
four  nickel  atoms  at  the  nearest  neighbor  distance  and  there  exists  an  aluminum  atom  at  site  A to 
lower  the  energy  of  site  A’.  Along  the  other  edge,  however,  both  B and  B’  sites  have  to  be 
occupied  by  two  aluminum  atoms  simultaneously  to  form  a unit  cell  and  such  an  occurrence  has 
a much  lower  probability  than  that  of  occupying  sites  A and  A’  in  sequence.  In  addition,  site  B’ 
has  only  two  Ni  atoms  at  the  nearest  neighbor  distance  and  the  probability  of  occupation  is  even 
lower  for  this  site.  Hence,  the  matrix  aluminum  atoms  tend  to  occupy  the  sites  immediately  next 
to  the  four-atom  edges.  This  argument  is  based  on  a two-dimensional  picture  but  it  is  trivial  to 
extend  the  logic  to  a three-dimensional  situation. 

The  aluminum  atoms  executing  random  motion  in  the  interface  region  arrive  at  sites  A and 
A’.  Once  the  aluminum  atoms  occupy  sites  A and  A’,  they  attract  other  aluminum  atoms  causing 
matrix  atom-depleted  regions  near  the  other  edges  of  dispersoid  boundary.  The  validity  of  this 
argument  is  supported  by  the  absence  of  such  unit  cell  formation  activity  by  the  matrix 
aluminum  atoms  when  the  symmetry  of  the  system  is  reduced  as  shown  in  the  results  of  Fig.  4 
obtained  with  the  dispersoid  lattice  rotated  by  15°  about  the  z-axis.  Nevertheless,  the  strongest 
evidence  for  the  above  argument  comes  from  the  results  of  45  ° rotation  shown  in  Fig.  5.  Here, 
the  five-atom  edges  are  compatible  with  matrix  lattice  and  the  unit  cell  formation  along  these 
edges  is  clearly  seen.  Such  site  preference  vanishes  at  elevated  temperature  as  shown  in  the 
simulation  results  of  Fig.  6 for  the  NijAl/Al  system  at  900K  with  the  dispersoid  rotated  by  30° 
where  a significant  reduction  of  the  matrix  atom  depletion  is  seen. 
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Fig.  7 shows  the  results  of  the  AljNb/AI  system  at  300K  obtained  without  rotating  the 
dispersoid.  The  system  under  this  situation  forms  a perfect  coherent  interface  without  any 
noticeable  lattice  strain  in  the  interface  region.  Fig.  8 shows  the  results  with  the  dispersoid 
rotated  by  30°  about  the  z-axis.  Here  again,  the  matrix  atoms  are  found  to  occupy  the  sites  to 
form  unit  cells  as  discussed  in  the  descriptions  of  the  model  shown  in  Fig.  3. 

The  most  unexpected  results  are  found  for  the  15°  and  45°  rotation  cases  as  seen  in  Figs. 

9 and  10.  The  two-dimensional  shape  of  the  dispersoid  in  these  two  figures  is  much  different 
from  that  appeared  in  the  results  of  Figs.  7 and  8.  Here,  the  dispersoid/matrix  bonding  is 


Figure  3.  Sites  favored  (A  &A’)  and 
less  favored  (B  & B’)  by  matrix 
atoms.  X isforNi  atz  = +a/2,  O 
for  Ni  at  z = 0 and  • for  A1  at  7 = 0. 
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Figure  4.  Simulation  results  of  the 
Ni3Al/Al  system  at  300K  with  the 
dispersoid  rotated  by  1 5°  about  the  z- 
axis. 
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Figure  5.  Simulation  results  of  the 
Ni3Al/AI  system  at  300K  with  the 
dispersoid  rotated  by  45°  about  the  z- 
axis 


Figure  6.  Simulation  results  of  the 
Ni3Al/Al  system  at  900K  with  the 
dispersoid  rotated  by  30°  about  the  z- 
axis 
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Figure  7,  Simulation  results  of  the 
A1  iNb/Al  system  at  300K  without 
rotating  the  dispersoid. 
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Figure  9,  Simulation  results  of  the 
AljNb/Al  composite  system  at  300K 
with  the  dispersoid  rotated  by  15° 
about  the  z-axis. 


Figure  8.  Simulation  results  of  the 
AhNb/Al  system  at  300K  with  the 
dispersoid  rotated  by  30°  about  the 
z-  axis. 
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Figure  10.  Simulation  results  of  the 
Al3Nb/Al  composite  system  at  3O0K 
with  the  dispersoid  rotated  by  45°  about 
the  z-axis. 


complete  as  if  they  were  interdiffused  into  each  other.  Such  a change  in  the  dispersoid  geometry 
may  or  may  not  be  related  to  the  fact  that  the  currently  assumed  LI  2 structure  of  A KNi  is  a 
metastable  phase.  No  explanation  is  available  at  the  moment  and  systematic  study  is  needed  to 
understand  the  results  of  these  figures.  Finally,  the  order  parameters  calculated  for  systems  with 
unrotated  dispersoids  over  a range  of  temperature  are  shown  in  Figs.  1 1 and  12  for  NhAl/Al  and 
AhNb/Al,  respectively.  It  is  seen  that  the  NijAl  dispersoid  is  thermally  stable  while  AKNb  is 
not.  This  is  puzzling  since  the  melting  temperature  is  considerably  higher  in  AljNb  than  NijAl 
(1953K  vs.  1673K).  The  melting  behavior  of  a particle  embedded  in  a matrix  can  be  strongly 
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Figure  11.  Temperature  dependence  of 
order  parameters  of  the  unrotated 
dispersoid  and  matrix  in  the  Ni3AI/AI 
system. 
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Figure  12.  Temperature  dependence  of 
order  parameters  of  the  unrotated 
dispersoid  and  matrix  in  the  AfvNbl/Al 
system. 


dependent  on  the  surrounding  environment  and  the  temperature  dependence  of  the  order 
parameter  needs  to  be  evaluated  under  various  situations  to  make  a general  statement. 

SUMMARY  AND  CONCLUSION 


The  simulation  results  of  the  present  study  show  the  general  trend  of  unit-cell  formation  in 
the  matrix  if  the  lattice  planes  of  the  matrix  and  dispersoid  are  appropriately  aligned.  The  case  of 
1 5°  and  45°  rotations  of  the  dispersiod  atoms  produce  matrix  atom-concentrated  and  depleted 
regions  near  the  dispersoid.  It  was  found  that  a perfectly  coherent  interface  forms  between 
Al3Nb  and  aluminum  matrix  under  appropriate  conditions.  Also  in  the  AhNb'Al  system,  the 
dispersoid  geometry  is  significantly  changed  while  its  wetting  with  the  matrix  remains  perfect 
for  particular  dispersoid  orientations.  The  underlying  mechanism  for  such  perfect  wetting  needs 
further  systematic  study  along  with  the  temperature  effects  and  stability  of  the  composite 
systems. 
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